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1. Introduction 

Most LGRBs are identified with massive stellar explosions, frequently in star-forming galax¬ 
ies. Many host galaxies (HGs) are observed and it is an open question yet whether the luminosity 
distribution of HGs, characterized by GRB events, is identical with those of the galaxies in general. 
Unfortunately, many HGs are fainter than the detection limit of optical observation. 

Nevertheless, we can give an upper limit of their luminosity [5]. Here the survival analysis can 
help to answer this question giving an unbiased estimate of the luminosity distribution. 

2. Mathematical summary 

The Survival analysis is a statistical procedure where some data in the sample of the mea¬ 
surements are censored (have only an upper/lower bound). Such data sets are also common in 
astrophysical problems. The detection limitations in the form of upper or lower limits are often 
called left or right censoring: the luminosity expressed in magnitude is a right censored event. The 
output of this procedure is an unbiased estimate of the survival function: 

S(t ) = Prob(T > t) 

where t is the measured value (magnitude in our cases) of the T random variable. The distribution 
function (Figure 1) which also called as lifetime distribution function is: 

F(t) = l-S{t) 

For the more astronomical context see Feigelson and Nelson [1]. 




Figure 1: Distributions calculated from the survival function with 95% confidence bounds in the I (left) and 
R (right) colors. 
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The basic Cox regression proportional hazards (Cox PH) model fits survival data depending 
on z co-variate to a hazard function (h(tlz)) in the form of 


^ = h(t\z) = h 0 -e^ 

where /J is an N-element unknown vector and ho(t) is called the baseline hazard, which is non- 
parametric. Primary interest lies in estimating the parameter using the partial likelihood function: 


N 


m=u 


rP'Zi 


i =o Lj€R(n) eP ' Zj 


where R(fi) is the risk set at t, [2], If the estimate of the baseline survival function So(t) is provided, 
then the estimate of the survival function for an individual with co-variates z^ may be obtained via: 


s{t\zk) = [ s 0 (or p {p ' Zk) 


3. Observational data 

For searching host galaxies more than 160 GRB positions were observed in the Preliminary 
Keck GRB Host Project. Most observations were carried out between 2005 and 2010 under a series 
of proposals focusing on host discovery and basic characterization (via the observers-frame optical 
color), and putting redshift constraints, in particular to rule out a large high-z fraction which was 
suggested in some early works [4]. We used the Preliminary Keck GRB Host Project Imaging 
Catalog with 164 host galaxies. The photometric systems were Johnson’s UBV and Cousins’ RI. 

The DEEP2 survey used the Keck II telescope to study the distant Universe. First the LRIS 
spectrograph was used in Phase 1 and a sample of 1000 galaxies studied within a limit of 1=24.5 
magnitude. Then in Phase 2 of the DEEP project was using the DEIMOS spectrograph to ob¬ 
tain spectra of 50,000 faint galaxies with redshift of z>0.7 and a limiting apparent magnitude of 
Rab=24. 1 mag. This spectrograph can observe in three bands (BRI). The survey covers an area 
of 2.8deg 2 divided into four separate fields observed [3], We assumed that all measured data are 
reliable namely not censoring data so we could use over than 49,000 galaxies in I band and about 
48,000 galaxies in R band. 

4. Survival analysis data 

We used the R and I magnitudes of Preliminary Keck GRB Host Project representing a suffi¬ 
ciently large homogeneous sample for survival analysis in R programming language with ’survival’ 
package [6]. We have 58 R and 38 I color observation having 14 and 7 lower limits (censored) for 
the magnitudes, respectively. We calculated both (R and I) survival functions from these magni¬ 
tudes. 

We used the DEEP2 galaxy catalog for making the comparison to the results of the survival 
analysis. The compared range was between about 21 and 24 magnitude in both colors because here 
are enough many data in both database. 

Testing the luminosity distributions no significant difference was found between HGs with PL 
or CPL LGRB spectra. The Cox regression was used to determine whether the properties of GRBs 
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(SWIFT BAT T90, Fluence, 1-sec Peak Photon Flux, Photon Index, XRT 11 Hour Flux, XRT 24 
Hour Flux and redshift) depend on the HG distribution. 

5. Results 

Figure 2 shows the difference between the cumulative distribution of the measured data and 
survival function in both color. It seems that the fainter are the objects more different are the two 
functions because the most censored data are among fainter points. 



Figure 2: The comparison of the distribution function obtained from the survival function (blue) with the 
cumulative measured data (red) of I color (left) and R color (right). 

Since the GRBs are connected to the standard candle luminosity function assumption [7, 8], 
we compared the HGs luminosity distribution with DEEP2 galaxies sample between about 21 and 
24 magnitude on both colors. Significant difference in the two cases were found. The probabilities 
of random differences are only 3 • 10~ 5 in I band and 9 • 10 5 in R band. The results are shown in 
Figure 3. We found that the probability is three times smaller in I color meaning that the GRB host 
galaxies are bluer than normal galaxies. 

Table 1 shows the result from the basic Cox regression. We didn’t find significant correlations 
between the GRB host brightness and some GRB physical parameters as the BAT T90, Fluence 
and Photon index of the GRBs. But we didn’t yet examine all off physical parameters of GRBs we 
would like to end this work. 

When we tested the CPL-PL groups difference in host galaxies with log-rank or Cochran- 
Mantel-Haenszel test we found there are no differences between the two classes, the probability 
being 0.73 in R and 0.5 in I band. 

Finally, almost the same results were found when we examined that if the redshifts are known, 
the probabilities are around 0.5. This result is very surprising because it has been expected that the 
redshifts of the fainter host galaxies would be less known. 
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Figure 3: The comparison of the distribution function obtained from the survival function I color and the 
DEEP2 I color distribution (left) and R color (right). The brightness distribution of the hosts are systemat¬ 
ically fainter than the DEEP2 galaxies, the probability of random difference is 9 ■ 10 s in R band and only 
3 • 10" 5 in I band. 


Variable 

I color 

R color 


i^cocff 

Prob 

Pcoeii 

Prob 

log(l T z) 

-1.125 

0.46 

-1.406 

0.21 

log(T90) 

0.182 

0.83 

-0.129 

0.85 

log (FI uence) 

-1.685 

0.17 

-0.540 

0.62 

log (PeakFlux) 

0.462 

0.68 

0.424 

0.63 

log(Photonindex) 

-2.675 

0.36 

1.923 

0.47 

log (XRT nh ) 

-1.621 

0.46 

1.072 

0.55 

log (XRT 24h ) 

1.578 

0.42 

-0.838 

0.60 


Table 1: Results of the Cox regressions. The Prob columns gives the probability, the dependence values are 
fully random. Seemingly, no remarkable correlations were obtained. 


6. Conclusion and Summary 

Comparing the cumulative distribution functions resulted in the survival analysis of the GRB 
hosts with that obtained from the DEEP2 sample in the R and I colors we recognized an offset of 
hosts towards fainter magnitudes. This offset is significant at the 9 • 10~ 5 level in the R and 3 • 10~ 5 
in the I color. Assuming the same spatial distribution of the GRB hosts and the DEEP2 sample 
galaxies this offset is explained by the lower absolute brightness and probably systematically lower 
mass of the GRB hosts. Moreover, the offset in the R seems to be smaller than in the I color 
indicating a bluer R-I color index of the GRB hosts. The bluer color index may indicate a higher 
star forming activity in the hosts then in the galaxies of the DEEP2 sample. 

Moreover, we studied dependence of the GRB host galaxies brightness distribution on the 
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GRB SWIFT data as SWIFT BAT T90, Fluence, 1-sec Peak Photon Flux, Photon Index, XRT 11 
Hour Flux, XRT 24 Hour Flux and redshift and it wasn’t found significant relationship between the 
data. Finally, no significant differences were found between HGs with PL or CPL GRB spectra or 
the redshift of host galaxies (redshift of GRBs) arc known or not. 
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